Comment on "Fluctuation-induced first-order transition p-wave superconductors" 
by Qi Li, D. Belitz, and J. Toner [Phys. Rev. B 79, 054514 (2009)] 
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In this Comment, we show that the paper by Qi Li, D. Belitz and J. Toner, published in Phys. Rev. 
B 79, 054514 (2009), contains an incomplete mean-field analysis of a simple model of Ginzburg- 
Landau type. The latter contains a stable non-unitary phase, which has is not been found in 
this study and is missing in the outlined picture of possible stable phases. Both mean-field and 
renormalization-group studies in this publication overlook previous papers, including those in Phys. 
Rev. B. As it turns, the authors re-derive and publish known results. 

PACS numbers: 64.60.ae,64.60.De, 74.20.De 



A paper by Qi Li et al. [1] is intended to present new 
results about global and local gauge effects on the phase 
transitions in unconventional (p-wave) superconductors. 
However, this article contains errors, misleading texts, 
and demonstrates a bold neglecting of previous papers @- 
5], where the same results have already been published. 
The study in Ref. [3is based on the same methods as 
those applied in Ref. [HI] (see also Ref. (fill): (i) mean- 
field (MF) approximation [fj , (ii) MF-like approximation 
[H, 0], and (hi) a renormalization group (RG) investiga- 
tion within the one-loop approximation combined with 
e-expansion to first order in e = (4 — d); d is the spa- 
tial dimensionality 0. As in previous studies [HU, the 
main tasks in Ref. if are MF derivation of the possi- 
ble phases and RG study of the so-called "fluctuation- 
induced weakly-first-order phase transition" Q in p-wave 
superconductors. 

Here we will present compelling arguments in favor of 
our point of view for essential disadvantages of the paper. 
Besides, we will correct some essential errors in Ref. [l[. 

The authors of Ref. [l[ perform a MF analysis of the 
free energy density 



f(i>)=t\^\ 2 +u\- 
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(cf. Eq. (3.1) in Ref. [l(), where t, u, and v are vertex 
(Landau) parameters, and ip — (ip\, V>2> V^) is a complex 
vector field (the order parameter of p-wave superconduc- 
tors). They use the suitable representation 



ij) = ipo (n cos 4> + ifa sin ( 



(2) 



by the modulus ipo = \'4>\, the auxiliary ("phase") angle 
(j), and the unit vectors n and fa (n.rh = cos 8). The 
unit vectors h and fa point out the directions of the real 
and imaginary vector parts of the field i/> = (ip + iif) ): 
n || ip' , and fa || if) . 



The authors of Ref. [l| correctly claim that f(if>) is 
minimized by h — fa and ipo = —t/2u, provided v > 0, 
u > 0, and t < 0. The minimum is f(t, u) = —t 2 /4u and, 
as mentioned in Ref. [l|, describes the unitary phase - 
the phase, characterized by the property ip x xjj* = 0. 
According to these authors, in the case of v > 0, this is 
the only available stable order (ip ^ 0). 

This outline of stable phases in superconductors with 
v > and described by model ([IJ is completely wrong, as 
we shall see in the remainder of this Comment. First, the 
description of the unitary phase is incomplete. Second, 
there exists another phase - a non-unitary phase, which 
is stable when v > 0. 

The careful analysis of the equations of state df /dipo — 
df /d9 = df /d<j) = inevitably leads to two topologically 
different domains of the unitary phase (i/j x t/>* =0), as 
given by the actual minimizing condition n — ±to, i.e., 
6 = irl, where I = 0, ±1, . . . (in Ref. [l|, this condition is 
presented in a wrong way: n = fa, i.e., 6 — 2nl). The 
two domains of the same unitary phase are distinguished 
by the parallel and antiparallel mutual orientations of n 
and fa. 

The stability matrix S = {Sij} = (d 2 f /dfiid^j), 
where /i = fa — (ipo,9,(f)), describes the stability of the 
possible phases (equilibria) with respect to variations of 
the order parameters (ipo,9,(j)) around their equilibrium 
values. For the unitary phase (if) x ip* = 0), this matrix 
has the form 



5' 



-At, 
0. 
0. 



(vt 2 /2u 2 ) sin 2 
0, 



20, 



(3) 
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Using Eq. and the existence conditions (i < 0, u > 0), 
one obtains the stability conditions u > 0, v > 0, t > 0, 
as indicated in Ref. [l|. 

Besides, there are additional circumstances, which 
should be emphasized. It is readily seen from Eq. §3§ 
that the unitary phase is marginally stable with respect 
to ^-fluctuations for angles cf> — irk/2 (k = 0, ±1,...), 
whereas the same phase has a marginal stability to- 
wards (^-fluctuations for any 8 (in both cases, t < 0, 
u > 0, and v > 0). As shown in previous works, 
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see [I, Q and references therein, these marginal sta- 
bilities towards the angles 9 and 4> reflect the global 
gauge invariance of the model (JI]) and the phenomenon 
of global symmetry breaking in the six-dimensional space 
(5fti/>i , ^sipi , . . . , §hp3, 9^3) °f the order parameter vector 
i/>. Thus the zeros of the matrix elements S22 and S33, 
produced by some values of 9 and <fi, exhibit the fun- 
damental global symmetry of model (TTJ) rather than an 
uncertainty in the stability degree of the unitary phase 
(the same is valid for the non-unitary phases, discussed 
below) UM- 

Ref. |lj reports only for a particular type of non- 
unitary phase (if) x tp* ^ 0), where h J_ m, i.e., 9 = 
7r(7 + l/2); existence and stability conditions: (u + v) > 0, 
v < 0, and t < 0. The result for 9 follows as a direct so- 
lution of the equations of state. Obviously, this phase 
has the property if)' _!_ if)" . In Sec. I-III as well as any- 
where in the paper [l| , there are no references to previous 
papers with MF results for the phases of model (fTJ), al- 
though this model is often used (see , and references 
therein). Moreover, the authors do not make any stip- 
ulation for some purposive incompleteness of their con- 
sideration; neither in Sec. III. A nor elsewhere. Thus the 
reader remains with the wrong conclusion that the non- 
unitary phase of type if) J_ is the only non- unitary 
phase (if) x if)* ^ 0) described by the model (JXJ) . 

In fact, the careful MF analysis |2j shows that the 
model ([1]) contains a stable non-unitary phase in more. 
This second non-unitary phase (2ndNUP) is given by 



t 

"2u' 



-k. 



(4) 



(k = 0, ±1, . . . ), the free energy minimum / = — i 2 /4w 
- the same as for the unitary phase, and the stability 
matrix 



S 



-At, 0, 
0, 0, 



0, 0, (2vt 2 /u 2 ) sin 2 8 



(5) 



The 2ndNUP has the same modulus if)Q and energy / 
as the unitary phase. The difference with the unitary 
phase is in the restriction <f> = (n/2)k as well as in the 
equilibrium values of angle 9. The Eq. §5§ shows that 
the second non-unitary phase is stable for t < 0, u > 0, 
and v > - stability conditions, which are identical to 
those for the unitary phase. In the limit 9 — » irl the 
structure and the stability properties of 2ndNUP coincide 
with those of the unitary phase. Having in mind the 
same free energies of these phases, one may conclude that 
the unitary phase and the second non-unitary phase are 
identical for 6* — 5- 7rZ. These results show that the u — v 
phase diagram outlined in Fig. 1 of Ref. [l[ is wrong and 
should be corrected along our instructions. 

The authors of Ref. [1| have missed to find the 2nd- 
NUP, contained in model (fTJ). On the other hand, they 
have missed to take advantage from paper (and ref- 
erences therein), where this phase is described in MF 



(in slightly different notations; see Sec. II and especially, 
Sec. II.D.4 in Ref. 0). In Refs. the possible u-v 

phase diagrams of the free energy (TTJ) have been outlined 
in details as a particular case of a quite more general 
description and more realistic variants of such diagrams. 
Thus the statement in the end of the Introduction of ar- 
ticle [l[ that the authors "derive the MF phase diagram" , 
without making any stipulation for a lack of original con- 
tribution, is incorrect and unacceptable. 

In Sec. III.B of Ref. fl, the Eq. (3.5) for the MF- 
like renormalization of the free energy f(ip) contains an 
essential error. In this equation, the parameter w should 
be related to the charge q by w ~ q 3 rather than by 
the wrong relationship w ~ q, given on the line below 
the equation. In order to persuade himself in this, one 
should have a look on the original article [7( or some other 
sources, for example, Ref. [6j. 

Besides, within the MF-like approximation, in which 
the spatial fluctuations of the field if) are neglected but 
the fluctuations of the magnetic induction in the entire 
Ginzburg-Landau fluctuation Hamiltonian are preserved 
, the parameters t and u acquire g-renormalizations, 
too. The latter are important for the description of es- 
sential properties of the weakly-hrst-order phase transi- 
tion Q . Note, that in accord with Eq. (2e) in Ref. Q , the 
parameter v does not possess q-renormalizations within 
the model ([T]) , neither in the one- loop RG treatment nor 
within the MF-like approximation. 

In Sec. III.C-D of Ref. fl, for the aim of study of 
spatial fluctuations of ip, the approximation of uniform 
field [^(x) w if)] is abandoned and the x- variations of 
ip are taken into account. Within this extended scheme, 
the critical fluctuations of the field V>(x) can be stud- 
ied with the help of RG. Here we briefly discuss the 
RG analysis of the "action" (2.3), performed in Ref. [l[ 
(Sec. III. C) within the one-loop approximation and first- 
order e-expansion. Note, that the mentioned level of ac- 
curacy in the application of RG is enough to reveal the 
main features of the fluctuation effects on the properties 
of the phase transition from the normal metal state to 
the p-wave superconducting state at the phase transition 
temperature, where the mean magnetic induction is equal 
to zero, and a particular type of fluctuation-induced 
weakly- first-order phase transition may appear [(| 0] • 

First, we note, that Eq. (2.3) does not present an ac- 
tion, as wrongly noted in Ref. [f|. Rather, Eq. (2.3) 
presents a "generalized free energy (functional)" (alias, a 
"fluctuation Hamiltonian" rather than an "action" ) . The 
same is valid for the quantity Sa, given by Eq. (4a) and 
describing the time-independent spatial fluctuations of 
the magnetic induction for the particular case of uniform 
if). 

Second, the RG analysis within the e-expansion per- 
formed in Ref. [fj is an entire repetition of the RG study 
in Ref. Q. Apart of some slight differences in the nota- 
tions, the RG Eqs. (3.7a)-(3.7c) are the same as the first 
three equations (3) in Ref. [3| . Note, that owing to a sim- 
ple misprint in Ref. 0, the number factor 12 in the last 
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term of the equation for du/dl should be 96. This mis- 
print is without any consequence for the results derived 
in Ref. A second note is needed to point out that in 
Ref. [|| we use known RG equations for the parameters 
c, q, and \x (Refs. @, H|)- These equations are exactly 
the same as Eqs. (3.7d)-(3.7f) in Ref. fl. They do not 
contain the parameter v and are known from the original 
paper [7j as well as from numerous further publications; 
see, e.g., Ref. jy, Q. 

Thus our RG equations, presented in Ref. Q and the 
RG equations, re-derived in Ref. [l[ are the same. Then 
it is not strange that the authors of Ref. [l[ have rediscov- 
ered our RG results about the possible fixed points and 
their stability properties, as well as about the possible 
types of critical behavior and first-order phase transition. 

For example, using Eqs. (3.11a) and (3.11b), the au- 
thors of Ref. [l| rediscover the couple of fixed points of 
type u± ^ 0, u± 0, v± — (notations of Ref. [3|, where 
u corresponds to v in Ref. [l[ ; the fixed point coordinates 
v± correspond to the crystal anisotropy parameter con- 
sidered there but neglected in Ref. [l[). Moreover, the 
condition for the existence of this couple of fixed points, 
m > mb = rib/ 2 = 211, derived for the first time in 
Ref. [H has also been re-derived in Ref. [l[ in a slightly 
different notation: 2m — n > n^!~ 420.9 (in a minor dif- 
ference with the result in Ref. [3j). Let us remind, that 
within the present problem, the symmetry index n is an 
integer equal to the total number of component of the 
complex vector ip, while m = n/2 is equal to the number 
of components of the real vectors tp' and i//'. Following 
Ref. [ij, in our discussion of MF results for the phases, 
we have used n — 2m — 6. But in the RG studies large 
integer n are considered, too. 

The main conclusions of Refs. [l[ and Q are the same: 
The critical behavior appears for n > n&. For n < rib, 
the phase transition from normal metal state to p-wave 
superconductivity is a fluctuation-induced weakly-first- 
order phase transition. Under suitable physical condi- 
tions [6J] , this weakly- first-order phase transition may ap- 



pear in the known superconductors, where n does not 
exceed 6. 

The difference in the e-analysis performed in these two 
papers is in the wider scope of Ref. [|[ , where the crystal 
anisotropy has been taken into account, and the effect of 
this anisotropy on the stability of the critical behavior 
has been investigated in details. The crystal anisotropy 
enhances the effect of first-order transition. For exam- 
ple, in the present of crystal anisotropy, the critical be- 
havior appear for m = 2n > 5494 (see Ref.Q). Cer- 
tainly, in real metals and metallic compounds, the crystal 
anisotropy effect on the p-wave superconductivity cannot 
be neglected. 

Moreover, the RG analysis in Ref. [l[ is a limited case of 
the RG study presented in Ref. Q • In Ref. [3| , the simul- 
taneous effect of crystal anisotropy and quenched impu- 
rities has been obtained by a generalization of the simple 
fluctuation model (2.1)-(2.3), given in Ref. [l|. The sim- 
ple check of article contents undoubtedly demonstrates 
that the RG equations in Ref. [l| are a limited case of 
those in Ref. Q and describe a simplified picture, where 
both crystal anisotropy and quenched impurity effects 
are ignored. Hence.the RG analysis and the respective 
RG results in Ref. [l| exactly follow from the RG study 
in Ref. [3] when the parameters, describing the crystal 
anisotropy and quenched disorder are neglected. 

Finally, we summarize our main findings. In Ref. [lj, 
the MF analysis is incomplete and a stable non-unitary 
phase is missed. In result, the domain of stability of this 
phase should be indicated in the u — v phase diagram, 
shown in Fig. (1). Excepting this error, the MF anal- 
ysis in Ref. [l[ repeats previous studies and re-discover 
known results (see Ref. Q and references therein). The 
parameter w in Sec. III. B should be corrected, as shown 
in this Comment. The RG study, performed in Ref. [l[ 
within the e expansion, totally repeats essential parts of 
the previous RG studies (Refs. % 0]) and, hence, re- 
derive known results. 
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